Voltage-probe and imaginary potential models for dephasing in a chaotic quantum dot 
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We compare two widely used models for dephasing in a chaotic quantum dot: The introduction 
of a fictitious voltage probe into the scattering matrix and the addition of an imaginary potential 
to the Hamiltonian. We identify the limit in which the two models are equivalent and compute the 
distribution of the conductance in that limit. Our analysis explains why previous treatments of 
dephasing gave different results. The distribution remains non-Gaussian for strong dephasing if the 
coupling of the quantum dot to the electron reservoirs is via ballistic single-mode point contacts, 
but becomes Gaussian if the coupling is via tunneling contacts. 
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I. INTRODUCTION 

Extensive theoretical work has provided a detailed 
description of the universal features of phase-coherent 
transport in classically chaotic systems, such as universal 
conductance fluctuations, weak localization, and a non- 
Gaussian conductance distribution. 1-12 The advances of 
submicron technology in the past decade have made these 
manifestations of quantum chaos in electronic transport 
accessible to experiment. 13-20 Although experiments on 
semiconductor quantum dots confirm the qualitative pre- 
dictions of the phase-coherent theory, a quantitative com- 
parison requires that loss of phase coherence be included 
into the theory. Two methods have been used for this 
purpose. 

The first method, originating from Buttiker, 21 is to in- 
clude a fictitious voltage probe into the scattering matrix. 
The voltage probe breaks phase coherence by removing 
electrons from the phase-coherent motion in the quan- 
tum dot, and subsequently reinjecting them without any 
phase relationship. The conductance G^ of the voltage 
probe (in units of 2e 2 /h) is set by the mean level spac- 
ing A in the quantum dot and the dephasing time r^, 
according to G^ = 2irh/T<f,A. This method was used in 
Refs. 7, 8, 13, and 20. The second method is to include a 
(spatially uniform) imaginary potential in the Hamilto- 
nian, equal to —ih/2T^. This method was used in Refs. 
9 and 11. 

The two methods have given very different results for 
the distribution of the conductance G, in particular in 
the case that the current through the quantum dot flows 
through single-mode point contacts. While the distri- 
bution P(G) becomes a delta peak at the classical con- 
ductance for very strong dephasing (t,£ — * 0) in the 
voltage-probe model, P{G) peaks at zero conductance 
in the imaginary potential model. It is the purpose of 
the present paper to reconcile the two methods, and to 
compute the conductance distribution in the limit that 
the two methods are equivalent. 

The origin of the differences lies with certain short- 
comings of each model. On the one hand, the imaginary 
potential model does not conserve the number of elec- 



trons. We will show how to correct for this, thereby 
resolving an ambiguity in the formulation of the model 
noted by McCann and Lerner. 11 On the other hand, the 
voltage-probe model describes spatially localized instead 
of spatially uniform dephasing. This is perfectly rea- 
sonable for dephasing by a real voltage probe, but it is 
not satisfactory if one wants a fictitious voltage probe 
to serve as a model for dephasing by inelastic processes 
occurring uniformly in space. A second deficiency of the 
voltage-probe model is that inelastic scattering requires 
a continuous tuning parameter t,),, while the number of 
modes in the voltage probe can take on integer val- 
ues only. Although the introduction of a tunnel barrier 
(transparency T^) in the voltage probe allows the conduc- 
tance G<f, = N^T^ to interpolate between integer values, 
the presence of two model parameters creates an ambi- 
guity: The conductance distribution depends on and 
separately, and not just on the product N^T^ set by 
the dephasing time. 

In this paper we present a version of the voltage-probe 
model that does not suffer from this ambiguity and that 
can be applied to dephasing processes occurring uni- 
formly in space. This version is equivalent to a particle- 
conserving imaginary potential model. We show that the 
absorbing term in the Hamiltonian can be replaced by an 
absorbing lead (the voltage probe) in the limit N$ — > oo, 
r0 — * at fixed G^ — N^T^. This is the "locally weak 
absorption limit" of Zirnbauer. 2 Both shortcomings of 
the voltage-probe model are cured: The limit iV^, — » oo 
together with ergodicity ensures spatial uniformity of the 
dephasing, while the conductance G^ is the only variable 
left to parameterize the dephasing rate. 

The outline of the paper is as follows. In Sec. II 
we recall the voltage-probe model and derive the limit 
N$ — > oo, — > at fixed N^T^ from the particle- 
conserving imaginary potential model. We then calculate 
the effect of dephasing on the conductance distribution in 
the case of single- mode point contacts (Sec. III). The dis- 
tribution narrows around the classical series conductance 
of the two point contacts when the dimensionless dephas- 
ing rate 7 = 2irh/T c f > A becomes 3> 1, but not precisely in 
the way which was computed in Refs. 7 and 8. In Sec. IV 
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we briefly consider the case of multi-mode point contacts 
(number of modes » 1), which is less interesting. We 
conclude in Sec. V. 



II. TWO MODELS FOR DEPHASING 

The system under consideration is shown in Fig. 1. It 
consists of a chaotic cavity, coupled by two point contacts 
(with iVi and iV" 2 propagating modes at the Fermi energy 
E F ) to source and drain reservoirs at voltages V\ and 
V 2 . A current I = I\ = —I 2 flows from source to drain. 
In the voltage-probe model, 21 a fictitious third lead {N$ 
modes) connects the cavity to a reservoir at voltage V^. 
Particle conservation is enforced by adjusting in such 
a way that no current is drawn [I^ = 0). The third lead 
contains a tunnel barrier, with a transmission probability 
which we assume to be the same for each mode. The 
scattering matrix S has dimension M = N\ + N 2 + N<p 
and can be written as 



Sll S12 Si<f, 
S21 S 2 2 S 2l p 
S(j>l S02 Sffjfj) 



(1) 



in terms of iVj x Nj reflection and transmission matrices 



Sij . Application of the relations 22 
h 



SkiN k 



yields the (dimensionless) conductance G 
(h/2e 2 )I/(Vi - V 2 ), 
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-G 



G\&G 
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G(f>i + G 
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(2) 
(3) 

(4) 



Using unitarity of S we may eliminate the conductance 
coefficients Gki which involve the voltage probe, 

G = -Gi 2 + (G* 11 + Gi 2 ){G 22 + gig) ^ 



Gn + G12 + G 2 i + G 



'22 



The remaining conductance coefficients are constructed 
from the matrix 



S = 



Sll S12 
S21 S22 



(6) 



which formally represents the scattering matrix of an ab- 
sorbing system. The first term in Eq. (5) would be the 
conductance if the voltage probe would truly absorb the 
electrons which enter it. The second term accounts for 
the electrons that are reinjected from the phase-breaking 
reservoir, thereby ensuring particle conservation in the 
voltage-probe model. 

The imaginary potential model relates S to a Hamil- 
tonian H with a spatially uniform, negative imaginary 




FIG. 1: Chaotic cavity, connected to current source and drain 
reservoirs (1 and 2), and to a voltage probe (0). The voltage 
probe contains a tunnel barrier (dotted line). The voltage V4, 
is adjusted such that 7^ = 0. 



potential — i-fA/4ir. As used in Refs. 9 and 11, it retains 
only the first term in Eq. (5), and therefore does not con- 
serve particles. We correct this by including the second 
term. We will now show that this particle-conserving 
imaginary potential model is equivalent to the voltage- 
probe model in the limit — > 00, — > 0, iV^r^ = 7. 

Our equivalence proof is based on the general 
relationship 23,24 

S= 1 -2mW\E F - H + iTrWW^)- 1 W (7) 

between the N x TV scattering matrix S (N — N\ + N 2 ) 
and the N' x N' Hamiltonian H (the limit N' — > 00 is 
taken later on). The Hamiltonian contains an imaginary 
potential, H^ v = H^ u — iS^jA/An, with H a Hermitian 
matrix. For a chaotic cavity, H is taken from the Gaus- 
sian ensemble of random matrix theory. 25 The N' x N 
matrix W has elements 24,26 



irW, 



fin 



■% n N'A (2T; 



i-2r-Vi-r„). (8) 



Here r n is the transmission probability of mode n in the 
leads and the energy A is the mean level spacing of H . 
We embed W into an N' x N' matrix by the definition 



W, 



fin 



for N < n < N\ and define 



7rVF 2 „ = ^H> 2 „ + ^„ 7 A/47r. 



(9) 



Substitution into Eq. (7) shows that S is an N x N sub- 
matrix of an N' x N 1 unitary matrix 



5=1- 2mW^{E F -H + mWW^W. 



(10) 



We have neglected the difference between and 
for 1 < (i < N, which is allowed in the limit N' — > 00. 
The matrix S is the scattering matrix of a cavity with 
three leads: Two real leads with Ni, N 2 modes, plus 
a fictitious lead with N' — N modes. The transmission 
probability T n of a mode in the fictitious lead follows 
from Eqs. (8) and (9), 



r„ = 



{N'A + ^WW 



00, 



(11) 
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where we have used that 7rW 2 „ = 7A/4-7T for N < n < 
N'. 

We conclude that the particle-conserving imaginary 
potential model and the voltage-probe model are equiv- 
alent in the limit iV = N' - N -> oo, T = j/N' -> 0, 
iV r = 7 (1 - N/N') -» 7. 



on the manifold of unitary (unitary symmetric) matrices 
for /3 = 2 (1). 

We now specialize to the case of ideal single-mode point 
contacts, Ni = N2 — 1 and Ti = T2 = 1. We seek the 
distribution of the 2x2 submatrix S defined in Eq. (6). 
We start with the polar decomposition of S, 



III. SINGLE-MODE POINT CONTACTS 

The effect of quantum-interference on the conductance 
is maximal if the point contacts which couple the chaotic 
cavity to the source and drain reservoirs have only a 
single propagating mode at the Fermi level. Then the 
sample-to-sample fluctuations of the conductance are of 
the same size as the average conductance itself. One thus 
needs the entire conductance distribution to characterize 
an ensemble of quantum dots. (An ensemble may be gen- 
erated by small variations in shape or in Fermi energy.) 

In the absence of dephasing, the conductance distribu- 
tion P(G) is strongly non-Gaussian. 3-6 For ideal point 
contacts (transmission probabilities Ti = r 2 = 1), one 
finds 4 - 5 



P(G) = i/3G^- 2 )/ 2 . 



(12) 



The symmetry parameter (3 = 2 (1) in the presence (ab- 
sence) of a time-reversal-symmetry breaking magnetic 
field. For high tunnel barriers (ri, r 2 <C 1), P(G) is max- 
imal for G = 0, and drops off oc G~ 3 / 2 for G > r^. 3 - 6 
In this section, we compute the conductance distribu- 
tion in the presence of dephasing, using the voltage-probe 
model in the limit — > 00, — > at fixed N^T^, in 
which it is equivalent to the current-conserving imaginary 
potential model. We focus on the case of ideal point con- 
tacts, and discuss the effect of tunnel barriers briefly at 
the end of the section. 

The scattering matrix S is distributed according to the 
Poisson kernel, 26-29 



u 
v 



it \/T — 



(14) 



where u and u' (v and v') are 2x2 (N^ x A^) unitary 
matrices, and t is a N$ x 2 matrix with all elements equal 
to zero except t nn = a/^TT, n = 1,2. In the presence of 



time-reversal symmetry, u' 



and 



of the polar decomposition (14) we have 







y/T=% 



In terms 



(15) 



The two parameters T\ and T 2 govern the strength of 
the absorption by the voltage probe. For T\,T2 — » the 
matrix S is unitary and there is no absorption, whereas 
for Ti , T2 — ► 1 the matrix S vanishes and the absorption 
is complete. Substitution of the invariant measure 12 



dS = \T 1 -T 2 f{T 1 T 2 )^ N *- 
x dudu'dvdv'dT 1 dT 2 



2-/3)/2 



(16) 



and the polar decomposition (14) into the Poisson kernel 
(13), yields the distribution of S in the form 



(T 1 T 2 )(/3^-2-/3)/2 



(17) 



\&et(l~v'vT)\f iN 4>+ 2 +P 

r = ^(i-r )(i-ttt). (is) 



P(S) 



(0M+2-0)/2 



1 det(l - fffft) 

V |dct(l-SSt)|/3M+2-,3 



(13) 



Since Eq. (Ill) is independent of u and it', the matrices 
u and v! are uniformly distributed in the unitary group, 
and the distribution of S is completely determined by the 
joint distribution P(Ti,T 2 ) of the absorption probabili- 
ties Ti and T 2 . 

It remains to perform the integral over v and v' in Eq. 
(III). This is a non-trivial calculation, which we describe 
in the appendix. The final result in the limit N^, — > 00, 
N 2 < n < M). The measure dS is the invariant measure — > at fixed 7 = iV^r^ is 



where V is a normalization constant, M = N\ + N2 + 
Nrj, is the dimension of S, and S is a diagonal matrix 
with diagonal elements S nn = VI — r„. Here T n is the 
transmission probability of mode n (T n = for Ni + 



P(T!,T 2 ) 



^Tf 4 T 2 - 4 exp 



\Ti-T 2 \ [7 2 (2-2 e T + 7 + 7 eT) 



- 7 (Ti + T 2 )(6 - 6e 7 + 4 7 + 2 7 e 7 + 7 2 ) + TiT 2 (24 - 24e 7 + I87 + 6 7 e 7 + 6 7 2 + 7 3 )] (19) 
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for (3 = 1 (presence of time-reversal symmetry) , and 

P(T U T 2 ) = iTf^exp^T^+T-T 1 )]^-^) 2 

x [ 7 4 (1 - 2e 7 + e 27 - 7 2 e 7 ) - 7 3 (Tx + T 2 )(4 - 8e 7 + 4e 27 + 2 7 - 2 7 e 7 - 2 7 2 e 7 - 7 3 e 7 ) 

+ 7 2 (T 2 + T 2 )(2 - 4e 7 + 2e 27 + 4 7 - 4 7 e 7 + 7 2 + 7 2 e 7 - 7 3 e 7 ) 

+ 7 2 r 1 T 2 (20 - 40e 7 + 20e 27 + I67 - 167e 7 + 4 7 2 - 8 7 2 e 7 - 4 7 3 e 7 - 7 4 e 7 ) 

- 77iT 2 (Ti + T 2 )(12 - 24e 7 + 12e 27 + 24 7 - 24 7 e 7 + 12 7 2 + 2 7 3 - 2 7 3 e 7 - 7 4 e 7 ) 

+ T 2 T 2 2 (12 - 24e 7 + 12e 27 + 24 7 - 24 7 e 7 + 24 7 2 - 12 7 2 e 7 + 87 s + 4 7 3 e 7 + 7 4 - 2 7 4 e 7 )] (20) 



for (3 = 2 (absence of time-reversal symmetry) . 

To relate the conductance G to T\, T 2 , u, and u', we 
substitute the polar decomposition of S into Eq. (5), with 
the result 

2 

2 

+ (T 1 +T 2 y 1 J2 lu.flu'j^TiTj. (21) 

Eqs. (Ill) and (21), together with the uniform distribu- 
tion of the 2x2 matrices u, v! over the unitary group, 
fully determine the distribution P(G) of the conductance 
of a chaotic cavity with two ideal single-mode point con- 
tacts. We parameterize u, u' in Euler angles and obtain 
P(G) as a four-dimensional integral, which we evaluate 
numerically. The distribution is plotted in Fig. 2 (solid 
curves) for several values of the dimensionlcss dephasing 
rate 7 = 2ttU/t^A. For 7 » l, 30 the conductance distri- 
bution becomes peaked around the classical conductance 
G=l/2, 

P(G)^^-(l + \x\-6p lX )e-W if 7 »l, (22) 

where x — 2j(3(G — 1/2). Notice that the distribution 
remains non-Gaussian for all values of 7. The limiting 
distribution (22) is plotted in Fig. 3, for (3 = 1 and 2. 
The average and variance of the conductance are 

(G) = i-^ /3l7 - 1 +0( 7 - 2 ), (23) 
varG = i(l + 2^ 1 ) 7 - 2 +0( 7 - 3 ). (24) 

The effect of dephasing was previously studied in Rcfs. 
7 and 8 for the case = 1 of an ideal voltage probe 
(without a tunnel barrier). The corresponding results 
are also shown in Fig. 2 (dotted curves). We see that the 
limit N<f, — > 00, — > results in narrower distributions 
at the same value of 7 = N^T^. In particular, the tails 
G — > and G — > 1 are strongly suppressed even for the 
smallest 7, in contrast with the case of the ideal volt- 
age probe. The physical reason for the difference is that 
keeping small and setting equal to 1 corresponds to 



dephasing which is not fully uniform in phase space, and 
therefore not as effective as the limit N^, — * 00, — > 0. 
For large 7, the difference vanishes, and the distribution 
(22) is recovered for an ideal voltage probe as well. (The 
fact that the conductance fluctuations around G = 1/2 
are non-Gaussian was overlooked in Refs. 7 and 8.) 

We have shown in the previous section that the 
voltage-probe model in the limit N$ — > 00, — > is 
equivalent to the particle-conserving imaginary potential 
model. The requirement of particle conservation is essen- 
tial. This is illustrated in Fig. 4, where we compare our 
results with those obtained from the imaginary potential 
model without enforcing conservation of particles. [This 
model corresponds to setting G = — G12 in Eq. (5) and 
was first solved in Ref. 3.] For 7>1, the imaginary po- 
tential without particle conservation yields a distribution 
which is maximal at G = 0, instead of a strongly peaked 
distribution around G = 1/2 [cf. Eq. (22)]. 

The first two moments of the conductance can be com- 
puted analytically from Eqs. (Ill) and (21). The result- 
ing expressions (which are too lengthy to report here) are 
plotted in Fig. 5. The markers at integer values of 7 are 
the results of the ideal voltage- probe model of Refs. 7 and 
8, where r^, = 1 and 7 = N$ = 0, 1, 2, . . .. The remark- 
able result 8 that (G) is the same for 7 = and 7 = 1 
is special for dephasing by a single-mode voltage probe: 
The present model with spatially uniform dephasing has 
a strictly monotonic increase of (G) with 7 for (3=1. 

Sofar we have considered ideal point contacts. Non- 
ideal point contacts (i.e. point contacts with tunnel bar- 
riers) correspond to Fi,r 2 < 1 in the distribution (13) 
of S. This case can be mapped onto that of ideal point 
contacts by the parameterization 26-28 

S = R + T(l- S'R^S'T, (25) 

where R and T = iy 1 — R? are diagonal matrices. The 
only nonzero elements of R are Rn = \fl — T\ and 
R22 = \/l ~ r 2 . The distribution of S' is given by the 
Poisson kernel (13) with I\ = T 2 = 1. Physically, S' 
is the scattering matrix of the quantum dot without the 
tunnel barriers in the point contacts, while R (T) is the 
reflection (transmission) matrix of the tunnel barriers in 
the absence of the quantum dot. 26 We may restrict the 
parameterization (25) to the 2x2 submatrix S, 

S = R + fil-S'R^S'f, (26) 



■5 



cu 2 



6=1 




// 0.5/' X 















0.2 0.4 0.6 0.8 1 
G 




FIG. 2: Solid curves: Conductance distributions of a quantum 
dot with two ideal single-mode point contacts, computed from 
Eqs. (Ill) and (21) for dephasing rates 7 = 0, 0.5, 1, 2, and 5. 
The top panel is for zero magnetic field (/? = 1), the bottom 
panel for broken time-reversal symmetry (/3 = 2). The dotted 
curves are the results of Refs. 7 and 8 for the model of an ideal 
voltage probe (without a tunnel barrier), in which dephasing 
is not fully uniform in phase space. For 7 = the two models 
coincide. The value 7 = 0.5 is not accessible in the model of 
an ideal voltage probe (because 7 = N^T,/, can take on only 
integer values if = 1). 




FIG. 3: The limiting conductance distribution (22) for 7 S> 1 
(solid curves). A Gaussian distribution with the same mean 
and variance is shown for comparison (dotted curves). 
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FIG. 4: Solid curves: Same as in Fig. 2, bottom panel. Dot- 
ted curves: Results of the imaginary potential model without 
particle conservation. 



where the matrices S', R, and T are the upper- left 2x2 
submatrices of 5", R, and T, respectively. The matrix S' 
has the distribution given by Eqs. (Ill) and (III). The 
matrices R and T are fixed, so the distribution of S fol- 
lows directly from Eq. (26). 

For strong dephasing (7 3> Fi,^) we find that the 
conductance distribution becomes a Gaussian with mean 



and variance given by 

rir 2 



(G) = 



varG = 



ri + r 2 
4rfrl(r? 



2rfr|(4//?-r 1 -r 2 ) 

7(ri + r 2 )3 
-ri-Tiil-rfra) 



/37(ri + r 2 )3 



(27) 
(28) 



The average conductance (G) is the classical series con- 
ductance of the two point-contact conductances T\ and 
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FIG. 5: Variance of the conductance as a function of the 
dephasing rate 7, for (3 = 1 (solid curve) and (3 = 2 (dotted 
curve), computed from Eqs. (Ill) and (21). The crosses ((3 = 
1) and squares ((3 = 2) at integer 7 result from the model of 
Refs. 7 and 8 with the ideal voltage probe. The inset shows 
the average conductance for /3 = 1. (For (3 = 2 the average is 
trivially equal to 1/2 for all 7 in both models.) 



IV Fluctuations around the classical conductance are of 
order 7 -1 / 2 . For ideal point contacts (Fi,r2 — > 1) the 
variance (28) vanishes. The higher-order fluctuations are 
non-Gaussian, described by Eq. (22). 

Again our result is entirely different from that 
of the imaginary potential model without particle 
conservation, 3 ' 11 where P(G) becomes sharply peaked at 
G = when 7 » Fi,r 2 . We have verified that we re- 
cover the results of Ref. 3 from our Eqs. (Ill) and (21) 
if we retain only the first term in Eq. (5), i.e. if we set 
G = —Gi2- The results of Ref. 11 are recovered if wc 
symmetrize this term, i.e. if we set G = — (G12 + G 2 i)/2. 
(This is different from — G12 if (3 = 2 and 7 ^ 0.) Once 
particle conservation is enforced, the imaginary potential 
model leads unambiguously to Eq. (III). 



IV. MULTI-MODE POINT CONTACTS 

In this section we consider the case N\ , N2 3> 1 of a 
large number of modes in the two point contacts. The 
conductance distribution is then a Gaussian, hence it suf- 
fices to compute the first two moments of G. We first 
consider ideal point contacts (Ti = T2 = 1), and discuss 
the effect of tunnel barriers at the end. 

For Ni,N 2 3> 1 the integration over the scattering 
matrix S with the probability distribution (13) can be 
done using the diagrammatic technique of Ref. 31. The 



result for the average of the conductance coefficients Gij 
is 



(Gij) = 



N,N 3 



Aij — 



n + iv r 

(N + 7VV)3 



+ & (3, l Aij-, 



SijNi 



N + NaT, 



(29) 



-,(30) 



up to terms of order iV _1 . (We recall that N = N± + 
A 2 .) For the covariances cov (Gij, Gki) = (GijGki) — 
(Gij)(G k i) we find 



cov (Gij, Gki) = AikAji + Sp.iAuAjk 



+ 



2N l N J N k NiN tP (N 4> + AQT|(1 - I» 
(3(N + 7V r )6 



(31) 



In order to find the average and variance of the conduc- 
tance in the presence of dephasing, we substitute Eqs. 
(IV) and (31) into Eq. (5). The result is 



(G) 



varG = 



A^ 2 



1 



3/31 



N \ N ■ 
2iV 1 2 iVf 
/3iV 2 (iV + 7) 2 ' 



(32) 
(33) 



with 7 = 7v r . 

Eq. (32) was previously obtained by Alciner and 
Larkin. 10 Eq. (33) for varG agrees with the interpolation 
formula of Baranger and Mello 7 . The present derivation 
shows that this interpolation formula is in fact a rigorous 
result of perturbation theory. [However, the interpola- 
tion formula of Ref. 7 for (G) differs from Eq. (32).] In 
the final expression for (G) and varG only the product 
N^T^ appears, although the moments of the conductance 
coefficients Gij depend on and separately. Appar- 
ently, in large- N perturbation theory the precise choice 
of N,p and in the voltage-probe model is irrelevant, the 
conductance distribution being determined by the prod- 
uct N^T^ only. For small dephasing rates 7 <C N, Eq. 
(IV) agrees with Efetov's result, 9 who used the imaginary 
potential model without enforcing particle conservation. 
However, for 7 > N, our result differs from that of Ref. 
9, indicating the importance of particle conservation once 
the dephasing rate 7 and the dimensionless escape rate 
N through the point contacts become comparable. 

We have carried out the same calculation for the case 
of non-ideal point contacts. The transmission probability 
of mode n is denoted by T n (n — 1, . . . , N\ corresponding 
to the first point contact, n = Ni + 1, . . . , N\ + N 2 to the 
second point contact). The result is 



(G) 
varG 



9i9i 



3/31- 



929? 



9 2 i92 



(34) 



g 2 (g + i) 

^9i92-9i9' 3 +929' 1 4 -93g' 1 4 ) 



+ 



pg 2 (g + j) 2 
6(g 2 1 g , 2 + 929'i 2 ) 2 



+ 



/3.9 4 (.9 + 7) 2 



Pg 4 (g + j) 

7igi(ggg£ + gggj) 
P9 3 (9 + 7) 2 ' 



7 



(35) 



iVi N t +N 2 

g p = J^K, g' p = r - 9 = 9i+9v (36) 



n=l 



n=l+Afi 



One can check that Eq. (IV) reduces to Eq. (IV) for ideal 
point contacts (when g p = N\, g' p = N 2 ). As in the case 
of single-mode point contacts, varG oc 7~ 2 for 7 » 1 
without tunnel barriers, while varG oc 7 _1 otherwise. 



V. CONCLUSION 

In summary, we have demonstrated the equivalence of 
two models for dephasing, the voltage-probe model and 
the imaginary potential model. In doing so we have cor- 
rected a number of shortcomings of each model, notably 
the non-uniformity of the dephasing in the voltage-probe 
model of Refs. 7 and 8 and the lack of particle conserva- 
tion in the imaginary potential model of Refs. 9 and 11. 
We have calculated the distribution of the conductance 
and shown that it peaks at the classical conductance for 
strong dephasing once particle conservation is enforced, 
thereby reconciling the contradictory results of Refs. 7 
and 8, on the one hand, and Refs. 9 and 11, on the other 
hand. We find that for ideal single-mode point contacts 
(no tunnel barriers), conductance fluctuations are non- 
Gaussian and oc T0 for strong dephasing (r^ — > 0). In 
the case of non-ideal point contacts (with tunnel barri- 
ers), fluctuations are larger (oc ^/t%) and Gaussian for 

T0 -> . 

The effect of dephasing becomes appreciable when the 
dimensionless dephasing rate 7 = 27t?i/t0A is of the 
same order as the dimensionless escape rate g = ^2 n T n 
through the two point contacts. For 7^3, the weak- 
localization correction SG = (G)((3 = 2) - (G)((3 = 1) 
and the conductance fluctuations are given by 30 



6G = aig/7 
varG = big/ j- 



b 2 (gh) 2 



0{ghf 



(37) 
(38) 



where ai, b\, and b 2 are numerical coefficients determined 
by Eqs. (Ill), (III), (IV), and (IV). For the special case 
of two single- mode point contacts, we have 



ai 



h = 



4r 2 r 2 



(ri + r 2 ) 4 

AT 2 T 2 2 {T\ - 



r 1 r2-r2r 2 ) 



/3(ri + r 2 ) 4 



(39) 



(40) 



The coefficient b 2 is only relevant if ri,T 2 k, 1, when 
h w (2 - Ti - r 2 )/4/3 < 1 and b 2 w (1 + 2 < 5 (31 )/16. At 
finite temperatures, in addition to dephasing, the effect 
of thermal smearing becomes important. 9 Since thermal 
smearing has no effect on the average conductance, the 
weak-localization correction SG provides an unambigu- 
ous way to find the dephasing rate 7. 

The fact that dephasing was not entirely uniform in 
phase space in the model of Refs. 7 and 8 leads to small 



but noticeable differences with the completely uniform 
description used here, in particular for the case of single- 
mode point contacts. The differences may result in a 
discrepancy A7 m 1 in the estimated value of the dimen- 
sionless dephasing rate 7, if the ideal voltage-probe model 
of Refs. 7 and 8 is used instead of the model presented 
here. A difference A7 s=a 1 is relevant, as experiments on 
semiconductor quantum dots can have dephasing rates 
as low as 7 w 2. 32 

Both the voltage-probe model and the imaginary po- 
tential model only provide an effective description of de- 
phasing. They cannot compete with a microscopic the- 
ory of inelastic scattering in quantum dots (see e.g. Refs. 
33 and 34). At this time, a microscopic theory for the 
effect of inelastic scattering on the conductance distribu- 
tion does not yet exist. For the time being, the model 
presented here may well be the most realistic description 
available. 
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CALCULATION OF P(T 1 ,T 2 ) 

We start the calculation of P(Ti,T 2 ) from the integral 
expression (III), in which we may replace the double in- 
tegral of v and v' by a single integral of the matrix v'v 
over the unitary group (for (3 = 2) or over the manifold 
of unitary symmetric matrices (for (3=1). We make a 
substitution of variables v'v — > w via 



v v 



\J\ — t 2 w(l — tw) 1 \]\ — t 2 . 



(41) 



The matrix r was defined in Eq. (18). One verifies that 
the matrix w is unitary (unitary symmetric for /3 = 1). 
The Jacobian of this transformation is 26-28 

,'dv'v\ V \dct(l-v'vT)\t 3N *+ 2 -t 3 , s 
det ( — ] = - — , (42) 



dw 



V dct(l-T 2 )(' 3Ar *+ 2 -/ 3 )/ 2 ' 



where V and V are normalization constants. This 
change of variables is a key step in the calculation, since 
the Jacobian (42) cancels the denominator of the inte- 
grand of Eq. (17) almost completely, 

P(Ti,T 2 ) = lj I ' d W T^-V\ Tl -T 2 f 

x [] (i + r,T J - 1 -r,)-(' 3Ar *+ 6 -' 3 )/ 2 

j=l,2 

x \\ T r 2t) - 2 1 det(l - rw)\ 2t) . (43) 



.7=1,2 
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We now consider separately the integral 

Ip = J dw|dct(l-™)| 2/3 

= J dw det(l — y/rwy/r) 13 det(l — y/rw^ 1 y/r). 

(44) 

Here we have used that r is a positive diagonal matrix. 
We now change variables ^/tw^^/t — > w^ 1 . If the ma- 
trix r were unitary, we could write 

Ip = j dw det(l-rwr) /3 det(l-w- 1 ) /3 , (45) 

in view of the invariance of the measure dw — dw. How- 
ever, t is not unitary. A theorem due to Weyl allows us 
to continue Eq. (45) analytically to arbitrary t. 35 

To evaluate Ip, we decompose w in eigenvectors and 
eigenphases, w — Ue l ®U\ where U is an orthogonal (uni- 
tary) matrix for (3 = 1 (2), and 6y = Sij6j, < 9j < 2n. 
The invariant measure dw reads 25 

dw = dU]J\e i6i -e i6 tfY[d6i. (46) 

i<j i 

After some algebraic manipulations, we arrive at 
r r N + 

ip = d/k... de N , n ^ - & f n (! - 

J J i<3 3 = 1 

1] [i-(i-r*)e w 'f J dUdetA?, (47) 



3=1 



where the 2x2 matrix ^4 is given by 



A 



4 - • — 

V l 3 



■ (48) 



The determinant of A is computed by a direct expan- 
sion. Since ^> 1, we may consider the matrix elements 
Uki as independent real (complex) Gaussian distributed 
variables with zero mean and variance l/N$ for (3 = 1 
(2). We write the result of the Gaussian integrations in 
terms of derivatives of a generating function Fp , 

[1 - (1 - T<t>Y 6i ] / dtf det A" - DpFp. (49) 

The generating function Fp depends on the variables x k , 
yk, and z k , where k = 1 for (3 = 1 and fc = 1, 2 for /? = 2, 

p p = n n^ 1 +** +/(^,y fc ,^)e^i, (50) 



/(a;,y,z) = (l + a ; + y)- 1 (l-r^ 



l + a;(l-2T 1 ) + y(l-2T 2 ) + zV7V^J • (51) 
The differential operator Dp reads 



£>2 



\^\d Xl +d Vl ) + N^ 2 d Zl d Zl , 



(52) 



a: 



^(^ Xl 9x 2 + d yi d V2 ) - i(9 Xl - <9<, 2 ) 2 



+ N^{-d Z2 d Z2 - -a zi a zi )(a xi + d yi ) 

+ N7 4 d Zl d Z2 d Z2 (3d Z2 ~2d Zl ). 



(53) 



The derivatives in Eq. (49) should be evaluated at x k — 
Vk = Zk = (k = 1,2). 

We are left with an integral over the phases Oj which 
is of the type 



I'p = f dBx...j d0 n JJ|e w '- 

i<j 



,ie j i/3 



n 







xl[(l-e-^fl[(a k -e^). (54) 

3=1 k=l 



The integrand is a product of secular determinants 
det(A — U) of a unitary matrix U. Integrals of this form 
were considered by Haake et al. 36 For j3 = 1 we can di- 
rectly apply the results in their paper, for (3 = 2 we need 
to extend their method to include a product of four sec- 
ular determinants. We find 



I'l 



(1 + n)(<+ 3 - 1) - (3 + n) Ql K +1 - 1) 

( fll -l) 3 (n-l) 
K +2 -l)(4 +2 -l) 
(d - l) 2 (a 2 - l) 2 

K+ 2 - fl2 1+2 )(n + 2) 
(at - l)(a 2 - 2)(ai - a 2 ) ' 



(55) 



(56) 



The desired integral Ip is obtained from Jg by substitu- 
tion of Eq. (V) with n = N^, a k = f(x k , yk, z k ) into Eqs. 
(V)-(V). Substitution of Ip into Eq. (43) then leads to 
the final result (III). 
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